Flow Ripple Reduction in Reciprocating Pumps by Multi-Phase Rectification

Reciprocating piezoelectric micropumps enable miniaturization in microfluidics for lab-on-a-chip applications such as organs-on-chips (OoC). However, achieving a steady flow when using these micropumps is a significant challenge because of flow ripples in the displaced liquid, especially at low frequencies or low flow rates (<50 µL/min). Although dampers are widely used for reducing ripples in a flow, their efficiency depends on the driving frequency of the pump. Here, we investigated multi-phase rectification as an approach to minimize ripples at low flow rates by connecting piezoelectric micropumps in parallel. The efficiency in ripple reduction was evaluated with an increasing number (n) of pumps connected in parallel, each actuated by an alternating voltage waveform with a phase difference of 2π/n (called multi-phase rectification) at a chosen frequency. We introduce a fluidic ripple factor (RFfl.), which is the ratio of the root mean square (RMS) value of the fluctuations present in the rectified output to the average fluctuation-free value of the discharge flow, as a metric to express the quality of the flow. The fluidic ripple factor was reduced by more than 90% by using three-phase rectification when compared to one-phase rectification in the 2–60 μL/min flow rate range. Analytical equations to estimate the fluidic ripple factor for a chosen number of pumps connected in parallel are presented, and we experimentally confirmed up to four pumps. The analysis shown can be used to design a frequency-independent multi-phase fluid rectifier to the desired ripple level in a flow for reciprocating pumps.


Introduction
There are mainly two ways of controlling fluid flow in microfluidics: active flow control and passive flow control. Active flow control requires an energy input from an actuator, and passive flow control does not require any actuators for fluid flow manipulation [1,2]. Active flow control is used in many applications, such as biomedical healthcare [3], food industry [4], inkjet printing [5], therapeutics [6], and soft robotics [7]. When compared to passive flow control, active flow control may provide greater control over the flow rate, a more precise flow, and closed-loop control [8][9][10]. Reciprocating pumps are one of the preferred pump types in active flow control for the delivery of gases and low-viscosity liquids (~1 mPa·s, e.g., water) [11,12]. They contain one or more pumping elements (pistons, plungers, or diaphragms) that reciprocate into and out of pumping chambers to produce the flow. The need for portable microfluidic platforms increased the interest in miniaturized reciprocating pumps, especially piezoelectric diaphragm pumps. These pumps are ideal for liquid flow rates of less than 50 µL/min and back pressures of up to 1.5 bar [13]. They do not require priming, have simple structures, and can be miniaturized and mass-produced [14].
Low-flow rates are necessary for applications such as organs-on-chips (OoC) technology, where organ functions are emulated with the aim of minimizing animal testing [15].
Low-flow rates are necessary for applications such as organs-on-chips (OoC) technology, where organ functions are emulated with the aim of minimizing animal testing [15]. The fluid flow through OoC environments needs to be stable to provide an appropriate fluid shear stress level on the cells in microfluidic chambers, e.g., below a threshold of 3 Pa [16], while mimicking vascular perfusion [17][18][19][20]. The purpose is to generate a steady flow and prevent adverse effects on cells during fluid flow, which could influence the stability and viability of the cells.
Piezoelectric diaphragm pumps are promising in terms of miniaturization [12]. The working principle and the schematic representation of the pump are shown in Figure 1. The pump consists of a chamber with a piezoelectrically actuated diaphragm and check valves at the inlet and the outlet (Figure 1a) to displace the fluid in one direction. In the suction mode, the diaphragm deforms upward to allow the fluid inside the chamber while the check-valve at the inlet is open and the check-valve at the outlet is closed. In the discharge mode, the diaphragm deforms downward to displace fluid out of the chamber, while the check valve at the inlet is closed and the check valve at the outlet is open. The functionality in fluidic circuit symbols is shown in Figure 1b. The flow generated by piezoelectric diaphragm pumps pulsates due to the reciprocating motion of the displacing diaphragm ( Figure 1c). The amplitude and frequency of the diaphragm are the main parameters for flow rate control. There are several methods reported in the literature to remove the pulsations in reciprocating pumps [21]. These methods include pump design modification by changing chamber outlet geometry [22,23], pulsation dampers, such as flexible tubing [24] and compliant membranes [25,26], controllable flow restriction by throttling with an orifice [27,28], and the multiplication of pistons [21,29,30]. Pump design modifications and dampers are often limited in their effectiveness and introduce compromises with regard to system size and response time [31]. In some cases, high pressures above 1 bar are required for dampers to be effective [32,33]. Furthermore, dampers cover only a narrow frequency range of a There are several methods reported in the literature to remove the pulsations in reciprocating pumps [21]. These methods include pump design modification by changing chamber outlet geometry [22,23], pulsation dampers, such as flexible tubing [24] and compliant membranes [25,26], controllable flow restriction by throttling with an orifice [27,28], and the multiplication of pistons [21,29,30]. Pump design modifications and dampers are often limited in their effectiveness and introduce compromises with regard to system size and response time [31]. In some cases, high pressures above 1 bar are required for dampers to be effective [32,33]. Furthermore, dampers cover only a narrow frequency range of a design and, thus, cannot reduce frequency ripples effectively over the entire flow rate range of the pump [34]. Controllable flow restriction is carried out by adding an extra proportional valve, which might produce intolerable high back pressures by increasing the resistance of the microfluidic system.
We propose that multi-phase rectification is a promising method to reduce the pulsations in reciprocating pumps. To the best of our knowledge, no research has been conducted to explore and quantify multi-phase rectification by introducing fluidic ripple factors, albeit there are a few two-phase rectification studies that reached improvements of only 35% [35] and 55% [36] for flow ripple reductions (having out-of-phase differences between pumping chambers).
In this study, we explore generalized multi-phase rectification schemes and their effect on flow rate and flow ripples. We give a theoretical analysis on multi-phase rectification; then, the obtained results were compared to experiments using off-the-shelf piezoelectric diaphragm pumps. The ripple factor value of the flow was used as a figure of merit for a steady flow with increasing numbers of parallel-phased micropumps. As a final step, we analyze and discuss the back pressure values and the influence of flow restriction.

Theory
This section gives analytical equations for the computation of fluid flow and flow ripples. The equations are based on the analogy of rectification in electrical circuits.

Rectification Analogy with Electrical Elements
The electrical current (I) flowing through an electrical resistor (R) for an applied voltage (V) is given by Ohm's law (V = I × R). Similarly, the fluid flow ( Q) through a resistive flow channel (R h ) for an applied pressure difference (∆P) between the inlet and the outlet is given by Hagen-Poiseuille's law (∆P = Q × R h ) [37]. This analogy is quite useful in designing purely resistive microfluidic circuits.
The process of converting alternating current (AC) to direct current (DC) is known as electrical rectification. Input current (I AC ) and average rectified output DC current (I DC ) over a time period of ω = 2π are given as [38] I AC = I max sin(ωt) (1) Since the π to 2π cycle of a 2π period of I DC is equal to zero for a half-wave rectifier, Equation (2) can be rewritten as Similarly, the root mean square value (I RMS ) is defined as The ripple factor (RF) is used to measure the ripple content in the rectified signal. The ripple factor is defined as the ratio of the root mean square (RMS) value of the AC component present in the rectified output to the average DC component of the rectified output signal when the load is purely resistive [39]: With the analogy of electrical current (I) in the electrical domain-to-flow rate (Q) in the fluidic domain, the following equation can be written for the flow rate in microfluidic rectifiers (with flow discharge only for 0 to π cycle of a 2π period in a single pump or the equivalent half-wave rectifier): where Q Stroke is the peak flow rate generated by one stroke of the actuator in the piezoelectric diaphragm pump. Here, valves and actuators are considered to be ideal for the sake of simplicity, i.e., valves have zero forward pressure drop, actuators have zero resistance, and no reverse flow or leakage occurs in the pumping system.

Multi-Phase Rectification
In multi-phase rectification, n micropumps are connected in parallel, as shown in Figure 2. Each individual pump is considered similar to the others in terms of the chamber, valve, and channel parameters, resulting in the same performance. The phase difference of actuation for each pump is 2π/n. A pictorial representation of the same is given in Table 1 with examples of up to five phase rectification. With the analogy of electrical current ( ) in the electrical domain-to-flow rate ( ) in the fluidic domain, the following equation can be written for the flow rate in microfluidic rectifiers (with flow discharge only for 0 to cycle of a 2 period in a single pump or the equivalent half-wave rectifier): where is the peak flow rate generated by one stroke of the actuator in the piezoelectric diaphragm pump. Here, valves and actuators are considered to be ideal for the sake of simplicity, i.e., valves have zero forward pressure drop, actuators have zero resistance, and no reverse flow or leakage occurs in the pumping system.

Multi-Phase Rectification
In multi-phase rectification, n micropumps are connected in parallel, as shown in Figure 2. Each individual pump is considered similar to the others in terms of the chamber, valve, and channel parameters, resulting in the same performance. The phase difference of actuation for each pump is 2 / . A pictorial representation of the same is given in Table 1 with examples of up to five phase rectification.   Table 1. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel for up to n = 5, five-phase rectification (flow profiles and descriptions from n = 6 (six-phase) to n = 9 (nine-phase) rectification are given in Appendix B). Each sinusoidal peak is the result of a piezoelectric stroke of each micropump (Figure 1). The different colors represent the flow profiles of different pumps. With an increase in the phase rectification (number of pumps), the strokes produced by pumps overlap, thus effectively increasing the Q DC with the number of pumps. The amount of overlap of the strokes also influences the Q DC−max − Q DC−min value: lower for an odd-phase rectifier between two adjacent, even-phase rectifiers.

Rectified Output Flow Profiles Description
Sensors 2023, 23, 6967. https://doi.org/10.3390/s23156967 www.mdpi.com/journal/s Table 1. The effect of rectification on the output flow rate for multi-phase rectification for p connected in parallel for up to = 5, five-phase rectification (flow profiles and descriptions fr = 6 (six-phase) to = 9 (nine-phase) rectification are given in Appendix B). Each sinusoidal p the result of a piezoelectric stroke of each micropump (Figure 1). The different colors represe flow profiles of different pumps. With an increase in the phase rectification (number of pump strokes produced by pumps overlap, thus effectively increasing the with the number of p The amount of overlap of the strokes also influences the − value: lower odd-phase rectifier between two adjacent, even-phase rectifiers.  Table 1. The effect of rectification on the output flow rate for multi-phase rectification for p connected in parallel for up to = 5, five-phase rectification (flow profiles and descriptions fr = 6 (six-phase) to = 9 (nine-phase) rectification are given in Appendix B). Each sinusoidal p the result of a piezoelectric stroke of each micropump (Figure 1). The different colors represe flow profiles of different pumps. With an increase in the phase rectification (number of pump strokes produced by pumps overlap, thus effectively increasing the with the number of p The amount of overlap of the strokes also influences the − value: lower odd-phase rectifier between two adjacent, even-phase rectifiers.  The effect of rectification on the output flow rate for multi-phase rectification for p connected in parallel for up to = 5, five-phase rectification (flow profiles and descriptions f = 6 (six-phase) to = 9 (nine-phase) rectification are given in Appendix B). Each sinusoidal p the result of a piezoelectric stroke of each micropump (Figure 1). The different colors represe flow profiles of different pumps. With an increase in the phase rectification (number of pump strokes produced by pumps overlap, thus effectively increasing the with the number of p The amount of overlap of the strokes also influences the − value: lower odd-phase rectifier between two adjacent, even-phase rectifiers. • Three pumps in parallel actuating at a 120 • phase difference. • Effective discharge flow rate has a cycloidal shape (dashed waveform), repeating itself at every 60 • phase difference, depending on the overlap of the individual fluid discharge profiles. Table 1. The effect of rectification on the output flow rate for multi-phase rectification for p connected in parallel for up to = 5, five-phase rectification (flow profiles and descriptions f = 6 (six-phase) to = 9 (nine-phase) rectification are given in Appendix B). Each sinusoidal p the result of a piezoelectric stroke of each micropump (Figure 1). The different colors represe flow profiles of different pumps. With an increase in the phase rectification (number of pump strokes produced by pumps overlap, thus effectively increasing the with the number of p The amount of overlap of the strokes also influences the − value: lower odd-phase rectifier between two adjacent, even-phase rectifiers.  For a period of T = 2π and an n-phase rectifier configuration, there will be n-pumps, with each pump being actuated at a phase difference (ϕ i ). Note that each pump discharges flow only in one half-cycle (00 to π) of the actuating full-cycle waveform:

Rectified Output Flow Profiles
The theoretically calculated values for Q DC , Q RMS , and RF f l. for up to nine-phase rectification with the pumps connected in parallel are shown in Figure 3. The Q DC and Q RMS values increase linearly with increasing numbers of phase-pump rectifiers because, with an increasing number of pumps, the number of pumps contributing to the effective output flow at a given phase also increases (see Table 1). The Q DC−max to Q DC−min values increase after every odd number of phase-pump rectifiers because an additional phase adds to the Q DC−min minimum value. The RF f l. value already drastically decreases by 96.6% for the three-phase rectifier, and the subsequent improvement is asymptotically low. For a nine-phase rectifier, the RF f l. value is 99.6% lower than that of the one-phase rectifier. Note that with the increase in the pumps, the amount of power needed adds up. A MATLAB code for calculating n pumps in parallel and a case study of up to three-phase rectification are given in Sections C and D, respectively. after every odd number of phase-pump rectifiers because an additional phase adds to the − minimum value. The . value already drastically decreases by 96.6% for the three-phase rectifier, and the subsequent improvement is asymptotically low. For a ninephase rectifier, the . value is 99.6% lower than that of the one-phase rectifier. Note that with the increase in the pumps, the amount of power needed adds up. A MATLAB code for calculating n pumps in parallel and a case study of up to three-phase rectification are given in Appendices C and D, respectively. . , dramatically decreases to 3.4% of that of the one-phase rectifier. The . value goes to 0.38% for a nine-phase rectifier. The . values for higher-phase rectifiers are given in Table 2 or can be calculated for any chosen number of phase rectifiers using the MATLAB code given in the supplementary information. ) for up to nine-phase rectifiers are shown for a chosen Q Stroke = 1 for every pump. The Q DC and Q RMS values increase as the phase rectification increases. For a three-phase rectifier, the ripple factor value, RF f l. , dramatically decreases to 3.4% of that of the one-phase rectifier. The RF f l. value goes to 0.38% for a nine-phase rectifier. The RF f l. values for higher-phase rectifiers are given in Table 2 or can be calculated for any chosen number of phase rectifiers using the MATLAB code given in the Supplementary Information.
The theoretically calculated flow rates and ripple factors for up to 100-phase rectifiers (using 100 pumps connected in parallel) are given in Table 2. It is important to highlight that the ripple factor of a rectifier having an odd number of pumps (after three pumps) is less than the rectifier with a neighboring even number of pumps, and it is the same as that of the rectifier with double the number of pumps. For example, the three-phase rectifier (n = 3) has a smaller ripple factor than the four-phase rectifier (n = 4), and it has the same ripple factor as the six-phase rectifier (n = 6). This relationship can be seen for each odd number of phases in multi-phase rectifiers. This is because of the way the fluid discharge profiles of each individual pump overlap with respect to the applied phase difference, as shown in Table 1 and Appendix B. The effective discharge profiles of each even-numbered phase rectifier (n = 2,4,6,. . .) have cycloidal repetitions equal to the applied phase difference. For odd-numbered rectifiers (n = 3,5,7,. . .), the cycloidal repetitions are equal to half of the applied phase differences. It is advantageous to operate odd numbers of micropumps in parallel while realizing multi-phase rectification.

Materials and Methods
The following components were used for testing: micropumps (mp6, Bartels Mikrotechnik GmbH, Dortmund, Germany), an arbitrary waveform generator in one-phase rectification case (RS PRO, RSDG 805, Northants, UK), and a microcontroller in a multi-phase rectification case to control the driving frequencies and phases (Nano Every, Arduino), with amplifiers having 100× gain (BD300, PiezoDrive, Shortland, Australia), a flow sensor for monitoring the flow rate data (MFS3, Elveflow, Paris, France), a microfluidic chip as a resistive load (Fluidic 155, microfluidic ChipShop GmbH, Jena, Germany), and a flow restrictor to test the low flow rates (1/4-28 Micro-Metering Valve, IDEX, Sulzberg, Germany). The mp6 micropumps have two piezoelectric chambers connected in series with the appropriate check valves typically operated 180 • out of phase. For our experiments, we used only one piezoelectric chamber and left the other chamber unconnected. The hydraulic resistance of the tubing and the flow rate sensor (≈1 × 10 3 Pa·s/m 3 and ≈3.4 × 10 10 Pa·s/m 3 , respectively) are considered negligible when compared to the microfluidic chip's hydraulic An experimental setup was prepared to test the multi-phase rectification schemes, as shown in Figure 4. The electrical control of the micropumps was realized by writing a sinusoidal signal generation code (with phase shifts) in Arduino IDE to create a frequency range of 0.25 to 10 Hz. In the one-phase rectification experiments, a frequency range of 0.1 Hz to 2 kHz was obtained by using the arbitrary waveform generator. The driving voltage was controlled by external potentiometers. Portable voltage amplifiers having two outputs were used for realizing high amplitudes for the piezoelectric micropump actuation. An experimental setup was prepared to test the multi-phase rectification schemes, as shown in Figure 4. The electrical control of the micropumps was realized by writing a sinusoidal signal generation code (with phase shifts) in Arduino IDE to create a frequency range of 0.25 to 10 Hz. In the one-phase rectification experiments, a frequency range of 0.1 Hz to 2 kHz was obtained by using the arbitrary waveform generator. The driving voltage was controlled by external potentiometers. Portable voltage amplifiers having two outputs were used for realizing high amplitudes for the piezoelectric micropump actuation. The inlet and outlet reservoirs were maintained at the same height for conducting the experiments. All pumps were kept at the same height as the inlet reservoir. The height of the outlet reservoir was raised only for the back-pressure experiments (the test setup is given in Appendix D). The sinusoidal signals and their phase shifts (PS) were generated in Arduino, shaped by a resistor-capacitor (RC) circuit, and amplified (Amp) 100× before connecting them to the micropumps. The setup was tested for one-phase, two-phase, three-phase, and four-phase rectifiers, respectively. The electrical connections of the components are shown in red, and the fluidic connections are shown in blue.

Performance of a One-Phase Rectifier
A one-phase rectifier was investigated by pumping deionized water. The study was carried out by activating only one pump with only one of the piezo chambers of an mp6 Bartels pump activated and leaving the other pump chamber unconnected, effectively using the pump as a one-phase rectifier. The one-phase rectifier was tested for flow over time at various frequencies and three different amplitudes (Figure 5a). Notably, the flow rate ripples are very high at low frequencies, with an increasing trend in the range of 0.1-5 Hz, followed by a decreasing trend between 5 Hz and 100 Hz. Beyond 100 Hz, low ripples were observed in the flow. Although an increase in the driving voltage increases the flow ripples, high-driving voltages are able to cover wider OoC flow rate ranges when . Schematic of an experimental setup for multi-phase flow rectification. The inlet and outlet reservoirs were maintained at the same height for conducting the experiments. All pumps were kept at the same height as the inlet reservoir. The height of the outlet reservoir was raised only for the back-pressure experiments (the test setup is given in Appendix D). The sinusoidal signals and their phase shifts (PS) were generated in Arduino, shaped by a resistor-capacitor (RC) circuit, and amplified (Amp) 100× before connecting them to the micropumps. The setup was tested for one-phase, two-phase, three-phase, and four-phase rectifiers, respectively. The electrical connections of the components are shown in red, and the fluidic connections are shown in blue.

Performance of a One-Phase Rectifier
A one-phase rectifier was investigated by pumping deionized water. The study was carried out by activating only one pump with only one of the piezo chambers of an mp6 Bartels pump activated and leaving the other pump chamber unconnected, effectively using the pump as a one-phase rectifier. The one-phase rectifier was tested for flow over time at various frequencies and three different amplitudes (Figure 5a). Notably, the flow rate ripples are very high at low frequencies, with an increasing trend in the range of 0.1-5 Hz, followed by a decreasing trend between 5 Hz and 100 Hz. Beyond 100 Hz, low ripples were observed in the flow. Although an increase in the driving voltage increases the flow ripples, high-driving voltages are able to cover wider OoC flow rate ranges when compared to low-driving voltages. Hence, the driving amplitude was fixed to 100 V to compare single-and multi-phase rectification by covering the 2-60 µL/min flow rate range when conducting the other experiments in this study. compared to low-driving voltages. Hence, the driving amplitude was fixed to 100 V to compare single-and multi-phase rectification by covering the 2-60 µL/min flow rate range when conducting the other experiments in this study.
(a) (b) Figure 5. The performance of a one-phase rectifier with increasing driving frequency for three different driving voltages is shown. (a) The flow rate values measured over time with increasing driving frequency for three different driving actuation amplitudes of a single-chamber mp6 micropump are plotted. The ripples are large at low frequencies when compared to high frequencies. Moreover, the acquired flow rates were also reduced by increasing the frequency. (b) The average flow rates and the corresponding electrical power consumption with the driving frequency of the pump for three operating voltages from the data of Figure 5a are plotted. The error bars represent the standard deviation from the average flow rate. Low ripple factor flow rates can also be obtained by increasing the frequency of the micropump, but electrical power consumption increases exponentially above 300 Hz for 100 V of driving amplitude. Note that frequency (x-axis) is shown on a logarithmic scale.
The average liquid flow rate derived from Figure 5a is plotted with frequency in Figure 5b, with the standard deviation shown as the error bars. The amount of power consumed for different frequencies is also plotted in the same graph. The average flow showed an asymmetrical top hat-like distribution with frequency. This means that low flow rates can also be obtained by increasing the frequency beyond (in this case) 100 Hz. However, the power consumption steeply increases with increasing frequency. Figure 5. The performance of a one-phase rectifier with increasing driving frequency for three different driving voltages is shown. (a) The flow rate values measured over time with increasing driving frequency for three different driving actuation amplitudes of a single-chamber mp6 micropump are plotted. The ripples are large at low frequencies when compared to high frequencies. Moreover, the acquired flow rates were also reduced by increasing the frequency. (b) The average flow rates and the corresponding electrical power consumption with the driving frequency of the pump for three operating voltages from the data of Figure 5a are plotted. The error bars represent the standard deviation from the average flow rate. Low ripple factor flow rates can also be obtained by increasing the frequency of the micropump, but electrical power consumption increases exponentially above 300 Hz for 100 V of driving amplitude. Note that frequency (x-axis) is shown on a logarithmic scale.
The average liquid flow rate derived from Figure 5a is plotted with frequency in Figure 5b, with the standard deviation shown as the error bars. The amount of power consumed for different frequencies is also plotted in the same graph. The average flow showed an asymmetrical top hat-like distribution with frequency. This means that low flow rates can also be obtained by increasing the frequency beyond (in this case) 100 Hz. However, the power consumption steeply increases with increasing frequency.
In dynamic applications of piezo actuators, the power consumption of the actuator increases linearly with the actuator's capacitance, driving voltage, and driving frequency. The following equations are used to calculate the power consumption of the piezo actuator [40,41]: where C is the capacitance of the actuator, V pp is the peak-peak driving voltage, and f is the driving frequency. The total current and total power consumption of the system can be calculated by using The individual pump's piezo actuator capacitance is defined by the manufacturer as 16 nF [42]. At 100 V pp and 1 Hz for driving voltage and frequency, an individual pump needs a 1.6 µA average current (≈3.6 µA I RMS ) for continuous operation with a sine wave. Increasing the driving frequency of the pump results in an increase in the flowing current (see Equation 13). For example, for a driving frequency of 100 Hz with the same driving voltage, the average current flow is 0.16 mA (≈0.36 mA I RMS ), and for 1 kHz, it is 1.6 mA (≈3.6 mA I RMS ). The calculated I RMS values are added to the 50 mA of current flowing through the amplifier, and then the I total values are obtained. The system power consumption is calculated from the voltage (24 V) and current (50 mA) values of the amplifier, resulting in 1.2 W until 300 Hz, with the increase in the pump's power consumption (with frequency) being negligible. However, the power consumption increased exponentially above 300 Hz (Figure 5b).
As seen in Figure 5, a lower frequency range (0.1-75 Hz) has high fluidic ripples when compared to a higher frequency range. Therefore, working at higher driving frequencies is advantageous to achieve low ripple factors, but at the cost of increased system power consumption [43] and decreased efficiency in terms of the check valves [44]. As a result, the flow rate starts to drop significantly.

Improvement in Ripple Factor
We show here that multi-phase rectification will allow low-power, ripple-free low flow rates. The results of the one-phase, two-phase, and three-phase rectifiers tested at an amplitude of 100 V and varying frequencies from 0.25 to 10 Hz are shown in Figure 6. The fluidic ripple factors and the percentage improvement of the ripples from one-phase rectifier are shown in Figure 7. The experimental results of the RF fl. values are plotted along with the theoretical values in Figure 8. An increase in the average flow rate can be seen with an increase in the number of phase rectifications, as predicted by the theory (Figure 3). By using three-phase rectification, i.e., three micropumps in parallel with 120 • phase shifts, an average improvement of 90.5% in ripple reduction was obtained for 0.25 to 10 Hz. This improvement is nearly independent of the driving frequency in the tested frequency range. When compared to the theoretical improvement of 96.6% in Figure 3 and Table 2, there is a~6% difference. The possible reason for this difference could be due to not having exactly identical pumps from the manufacturer. Each micropump has its own tolerance for the pumping membrane and valve thicknesses, thus affecting the stiffness of the pumps. In the theoretical study, the flow profiles are ideal sinusoids since the transfer function (H(s) = pressure/flow rate) is assumed as 1 when the flow profiles are obtained. On the other hand, we observe that the flow profiles are more like distorted sinusoids in the experiments. Flow profiles are affected by at least three parameters that are neglected in the theoretical study: the threshold pressure needed to open the check valves (Figure 1), fluid inertia, and the compliance of the silicone tubings. These parameters influence the shape of the flow profiles so as to create distorted sinusoidal waveforms when compared to the theoretically described sinusoidal waveform (Tables A3 and 3 in Appendix E). The threshold pressure required for opening the valves and fluid inertia effects are observed at the accelerating part of the flow profile, creating a delay when compared to the theoretically obtained flow profiles. Threshold pressure depends on the valve design [45,46], and it has a similar effect on a conducting fluid as threshold voltage does on diodes [47]. Fluid inertia influences the motion of the valves by adding mass to the system [48,49] and creates a delay at the accelerating part of the flow profile. The capacitive behavior of flexible tubing creates compliance in the microfluidic system, which has already been studied in microfluidics [24,50,51]. Although the capacitive effect has a major influence on flow profiles and flow ripple values, we observed that its influence on flow ripple improvements is limited to an average of~6% in the frequency range of 0.25 to 10 Hz. The delay in the descending part of the flow profile can be up to 1 s, affecting the time constant of the output flow profile (Table 3). Due to the capacitive effect, flow rates do not reach a zero level in multi-phase rectifiers, even at very low frequencies.
in the theoretical study: the threshold pressure needed to open the check valves ( Figure  1), fluid inertia, and the compliance of the silicone tubings. These parameters influence the shape of the flow profiles so as to create distorted sinusoidal waveforms when compared to the theoretically described sinusoidal waveform (Tables 3 and A3 in Appendix E). The threshold pressure required for opening the valves and fluid inertia effects are observed at the accelerating part of the flow profile, creating a delay when compared to the theoretically obtained flow profiles. Threshold pressure depends on the valve design [45,46], and it has a similar effect on a conducting fluid as threshold voltage does on diodes [47]. Fluid inertia influences the motion of the valves by adding mass to the system [48,49] and creates a delay at the accelerating part of the flow profile. The capacitive behavior of flexible tubing creates compliance in the microfluidic system, which has already been studied in microfluidics [24,50,51]. Although the capacitive effect has a major influence on flow profiles and flow ripple values, we observed that its influence on flow ripple improvements is limited to an average of ~6% in the frequency range of 0.25 to 10 Hz. The delay in the descending part of the flow profile can be up to 1 s, affecting the time constant of the output flow profile (Table 3). Due to the capacitive effect, flow rates do not reach a zero level in multi-phase rectifiers, even at very low frequencies. Figure 6. The influence of rectification on ripple reduction from one-phase to three-phases at a lower frequency range (0.25 to 10 Hz) is shown. The measurement durations are 20 s for each frequency. A one-phase rectifier uses one chamber of the mp6 pump, the two-phase rectifier uses a single chamber from two mp6 pumps connected in parallel, and similarly, a three-phase rectifier uses three mp6 pumps. Note that the flow rate also increases with the rectification level, as predicted and shown in Figure 3. Figure 6. The influence of rectification on ripple reduction from one-phase to three-phases at a lower frequency range (0.25 to 10 Hz) is shown. The measurement durations are 20 s for each frequency. A one-phase rectifier uses one chamber of the mp6 pump, the two-phase rectifier uses a single chamber from two mp6 pumps connected in parallel, and similarly, a three-phase rectifier uses three mp6 pumps. Note that the flow rate also increases with the rectification level, as predicted and shown in Figure 3.    Figure 7. The fluidic ripple factors and ripple factor improvement with increasing phase rectification are calculated from the data in Figure 6 and were analyzed. (a) Fluidic ripple factors of one-, two-, three-, and four-phase rectifiers, respectively, are shown. (b) The percentage improvement of the fluidic ripple factors ( .) of the two-phase and three-phase rectifiers compared to the one-phase rectifier is shown. Three-phase rectification has the lowest ripples and an average improvement of 90.5% when compared to one-phase rectification in the frequency range 0.25-10 Hz.   Table 3. Comparison of the theoretical and experimental flow profiles and fluidic ripple factors obtained for 0.25 Hz of actuation. The experimentally obtained flow profile shapes deviate from the theoretical sinusoidal waveform because of the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered in the theory. Therefore, the improvement in the experimental fluidic ripple factor also deviates from the expected theoretical improvement. See Table A3 in Appendix E for further analysis.  Effective discharge flow rate has a cycloidal shape (dashed waveform), repeating itself at every 36° phase difference. Table 3. Comparison of the theoretical and experimental flow profiles and fluidic ripple factors obtained for 0.25 Hz of actuation. The experimentally obtained flow profile shapes deviate from the theoretical sinusoidal waveform because of the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered in the theory. Therefore, the improvement in the experimental fluidic ripple factor also deviates from the expected theoretical improvement. See Table A3 in Appendix E for further analysis.

Estimated Flow Profiles by Theoretical Study Obtained Flow Profiles in the Experiments
Onephase Twophase Threephase 2 of 8  Five pumps in parallel actuating at a 72° phase difference.


Effective discharge flow rate has a cycloidal shape (dashed waveform), repeating itself at every 36° phase difference. Table 3. Comparison of the theoretical and experimental flow profiles and fluidic ripple factors obtained for 0.25 Hz of actuation. The experimentally obtained flow profile shapes deviate from the theoretical sinusoidal waveform because of the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered in the theory. Therefore, the improvement in the experimental fluidic ripple factor also deviates from the expected theoretical improvement. See Table A3 in Appendix E for further analysis.

Estimated Flow Profiles by Theoretical Study
Obtained Flow Profiles in the Experiments  Effective discharge flow rate has a cycloidal shape (dashed waveform), repeating itself at every 36° phase difference. Table 3. Comparison of the theoretical and experimental flow profiles and fluidic ripple factors obtained for 0.25 Hz of actuation. The experimentally obtained flow profile shapes deviate from the theoretical sinusoidal waveform because of the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered in the theory. Therefore, the improvement in the experimental fluidic ripple factor also deviates from the expected theoretical improvement. See Table A3 in Appendix E for further analysis.

Obtained Flow Profiles in the Experiments
Onephase Twophase Threephase 2 of 8  Five pumps in parallel actuating at a 72° phase difference.


Effective discharge flow rate has a cycloidal shape (dashed waveform), repeating itself at every 36° phase difference. Table 3. Comparison of the theoretical and experimental flow profiles and fluidic ripple factors obtained for 0.25 Hz of actuation. The experimentally obtained flow profile shapes deviate from the theoretical sinusoidal waveform because of the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered in the theory. Therefore, the improvement in the experimental fluidic ripple factor also deviates from the expected theoretical improvement. See Table A3 in Appendix E for further analysis.

Estimated Flow Profiles by Theoretical Study
Obtained Flow Profiles in the Experiments 0.485 60% 0.542 55%  Figures 5a and 6 show that the flow profiles of the one-phase rectifier return to zero below 1 Hz. This indicates a fluidic time constant of approximately 1 s in our experimental setup. The flow profiles reach zero for excitation frequencies lower than 1 Hz (1 s fluidic time constant). The flow profiles do not reach zero for excitation frequencies higher than 1 Hz. For very high excitation frequencies, the flow profiles do not reach zero; the fluidic system is too slow to catch up, and, hence, the flow rate is smoothed out. Therefore, we expect the phase rectification method to be very effective below 1 Hz (above the fluidic time constant of the system). Figures 6 and 7a,b clearly show that rectification is effective at least up to 10 Hz. This indicates rectification is effective until it matches the smoothness obtained by high excitation frequencies. The fluidic time constant of the system is due to compliant elements like pump membrane capacitance and silicone tubings. If non-compliant tubes are used, e.g., stainless steel tubes, the fluidic time constant will be shorter. The close-up flow profiles, the calculated fluidic ripple factors, and the flow oscillations are shown and compared with the theoretical flow profiles for 0.25 Hz in Table 3. These profiles are also merged and compared in Appendix E to analyze the deviations for each rectifier. Similar results for the 0.5 and 1 Hz frequencies are also shown in Sections F and G.

Influence of Multi-Phase Rectification on Back Pressure
There will be a back pressure exerted on the micropumps due to the flow resistance in the microfluidic circuit, thus influencing the pumping flow rate. Therefore, we studied the influence of back pressure on the micropump. In order to generate back pressure in the micropumps, the height of the outlet reservoir was raised above the inlet reservoir level (see Figure 4), thus increasing the hydrostatic back pressure. The average flow rate, Q, was measured as a function of the back pressure, P, for three types of rectifiers at 1 Hz, as per Figure 9. During the experiment, the microfluidic flow resistance was kept constant, and the height of the outlet reservoir (along with the microfluidic chip used as fluidic load resistance) was increased in steps of 100 mm to observe the effect of increasing hydrostatic back pressure on the average flow rate. The average flow rate linearly decreased with increasing back pressure, as expected. The one-phase and two-phase rectifiers showed similar behavior and could withstand a back pressure of about 40 mbar when the flow rate became zero. In contrast, the three-phase rectifier could withstand higher back pressure: up to 60 mbar. For a given flow rate, three-phase rectification can handle higher back pressure (P 3 >> P 1 & P 2 ), as shown in Figure 9. This is due to the distributed pumping load among three parallel pumps, as shown by the overlapping phases (Table 1). As the rectification level increases, the phase overlapping also increases, indicating the back pressure distributed over the pumps connected in parallel. Therefore, the back-pressure capability of the system increases with increasing fluidic rectification. The influence of increasing frequency on back pressure for rectification is shown in Appendix H. As the average flow rate of the micropump increases with frequency (Figure 6), the higher the micropump frequency, the greater the amount of back pressure the given rectifier can withstand. became zero. In contrast, the three-phase rectifier could withstand higher back pressure: up to 60 mbar. For a given flow rate, three-phase rectification can handle higher back pressure (P3 >> P1 & P2), as shown in Figure 9. This is due to the distributed pumping load among three parallel pumps, as shown by the overlapping phases (Table 1). As the rectification level increases, the phase overlapping also increases, indicating the back pressure distributed over the pumps connected in parallel. Therefore, the back-pressure capability of the system increases with increasing fluidic rectification. The influence of increasing frequency on back pressure for rectification is shown in Appendix H. As the average flow rate of the micropump increases with frequency ( Figure 6), the higher the micropump frequency, the greater the amount of back pressure the given rectifier can withstand. Figure 9. The average flow rate with back pressure for one-phase, two-phase, and three-phase rectifier mechanisms with micropump(s) actuated at 100 V and 1 Hz operating conditions. The dotted line indicates that for a given average flow rate, the three-phase rectifier can withstand higher back pressure (P3) compared to two-phase (P2) and one-phase (P1) rectifiers. The values of P1 and P2 are lower than P3 and are close together because there is no phase overlap for one-phase or two-phase rectifiers, whereas, for three-phase rectifiers, there is either a phase overlap indicating a back-pressure load distribution between pumps or a higher piezo stroke (see Table 1).

Obtaining Very Low Flow Rates
Changing the height of the reservoirs to obtain very low flow rates (Figure 9) is not a practical solution; therefore, a micro-metering valve (a flow restrictor) was used after the micropumps. The micro-metering valve limits the cross-sectional area of the flow channel, leading to an increase in the back pressure and a decrease in the flow rate. The height of Figure 9. The average flow rate with back pressure for one-phase, two-phase, and three-phase rectifier mechanisms with micropump(s) actuated at 100 V and 1 Hz operating conditions. The dotted line indicates that for a given average flow rate, the three-phase rectifier can withstand higher back pressure (P 3 ) compared to two-phase (P 2 ) and one-phase (P 1 ) rectifiers. The values of P 1 and P 2 are lower than P 3 and are close together because there is no phase overlap for one-phase or two-phase rectifiers, whereas, for three-phase rectifiers, there is either a phase overlap indicating a back-pressure load distribution between pumps or a higher piezo stroke (see Table 1).

Obtaining Very Low Flow Rates
Changing the height of the reservoirs to obtain very low flow rates (Figure 9) is not a practical solution; therefore, a micro-metering valve (a flow restrictor) was used after the micropumps. The micro-metering valve limits the cross-sectional area of the flow channel, leading to an increase in the back pressure and a decrease in the flow rate. The height of the inlet reservoir, micropumps, micro-metering valve, and outlet reservoir were kept at the same level. The micropump actuation was fixed at a driving voltage of 100 V and a frequency of 1 Hz. The flow rates (2, 5, and 15 µL/min) obtained by adjusting the micro-metering valve are shown in Figure 10a, and their corresponding ripple factors are shown in Figure 10b. At a flow rate of 15 µL/min, the ripple factor improvement for the three-phase rectifier is about 93% when compared to the one-phase rectifier. Although this improvement decreases to 82% and 50% at 5 and 2 µL/min flow rates, respectively, the reduction in flow ripples by the three-phase rectification is still significant. However, at very low flow rates, the micropump parts, such as the passive check valves, might become damaged due to increased back pressure. A bypass channel after the micropump connected to the inlet reservoir could minimize the damage to micropumps. the inlet reservoir, micropumps, micro-metering valve, and outlet reservoir were kept at the same level. The micropump actuation was fixed at a driving voltage of 100 V and a frequency of 1 Hz. The flow rates (2, 5, and 15 µL/min) obtained by adjusting the micrometering valve are shown in Figure 10a, and their corresponding ripple factors are shown in Figure 10b. At a flow rate of 15 µL/min, the ripple factor improvement for the threephase rectifier is about 93% when compared to the one-phase rectifier. Although this improvement decreases to 82% and 50% at 5 and 2 µL/min flow rates, respectively, the reduction in flow ripples by the three-phase rectification is still significant. However, at very low flow rates, the micropump parts, such as the passive check valves, might become damaged due to increased back pressure. A bypass channel after the micropump connected to the inlet reservoir could minimize the damage to micropumps.

Conclusions
The ripples in a flowing fluid generated by active reciprocating pumps, like piezoelectric micropumps, can be eliminated by multi-phase rectification. More than 90% of flow ripples are eliminated by connecting three pumps in parallel and operating them at a 120° phase difference. Connecting n pumps in parallel and actuating them at a phase difference of 360°/n will asymptotically minimize the ripple factors in the fluid flow. Experiments were performed using off-the-shelf commercial piezoelectric micropumps operating in one-phase, two-phase (180°), and three-phase (120°), and between an average flow rate of

Conclusions
The ripples in a flowing fluid generated by active reciprocating pumps, like piezoelectric micropumps, can be eliminated by multi-phase rectification. More than 90% of flow ripples are eliminated by connecting three pumps in parallel and operating them at a 120 • phase difference. Connecting n pumps in parallel and actuating them at a phase difference of 360 • /n will asymptotically minimize the ripple factors in the fluid flow. Experiments were performed using off-the-shelf commercial piezoelectric micropumps operating in one-phase, two-phase (180 • ), and three-phase (120 • ), and between an average flow rate of 2 µL/min and 60 µL/min. A single reciprocating pump actuated at higher frequencies could also eliminate ripples in the flow, but this results in a higher flow rate or higher power consumption. Ripples can be eliminated at a lower flow rate with a single pump by actuating at a low voltage and a high frequency. But this results in less control over the back pressure. Connecting pumps in parallel and operating them at a phase difference eliminates ripples independent of the pump actuating frequency, increases the effective flow rate, and enables the pump system to withstand higher back pressure. An odd number of pumps is more effective at eliminating ripples than an even number of parallel pumps. Every additional pump used consumes an extra voltage amplifier power. Therefore, a choice has to be made on how many micropumps to use based on the allowed flow ripple factor and the power consumption. If performance, cost, and complexity are to be balanced, then a three-phase rectifier is recommended. To conclude, multi-phase fluid rectification opens an opportunity of using reciprocating pumps for a ripple-free fluid flow, independent of actuating frequency, especially for the low flow rates (<50 µL/min) needed for many organs-on-chip applications.
Supplementary Materials: The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s23156967/s23156967/s1, A Matlab code to estimate the ripple factor for n-phase rectifier.

Data Availability Statement:
The data that support the findings of this study are available from the corresponding authors upon request. Appendix B Table A1. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from six to nine phases (continuation of Table 1).

Rectified utput Flow Profiles Description
• The system consists of six parallel pumps with a 60° phase difference.
• Discharge flow has a cycloidal shape, repeating itself for each 60° phase. Appendix B Table A1. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from six to nine phases (continuation of Table 1).

3, 23, 6967 3 of 8
Appendix B Table A1. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from six to nine phases (continuation of Table 1). Appendix B Table A1. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from six to nine phases (continuation of Table 1).

Rectified Output Flow Profiles
• Appendix B Table A1. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from six to nine phases (continuation of Table 1).

Rectified Output Flow Profiles
• The system consists of nine parallel pumps with a 40 • phase difference. • Discharge flow has a cycloidal shape, repeating itself for each 20 • phase. Table A2. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from one to nine phases next to each other (for ease of comparison).   Table A2. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from one to nine phases next to each other (for ease of comparison).   Table A2. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from one to nine phases next to each other (for ease of comparison).

Rectified Output Flow Profiles
Sensors 2023, 23 Table A2. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from one to nine phases next to each other (for ease of comparison).

Rectified Output Flow Profiles
Sensors 2023, 23 Table A2. The effect of rectification on the output flow rate for multi-phase rectification for pumps connected in parallel, from one to nine phases next to each other (for ease of comparison).

Rectified Output Flow Profiles
Sensors 2023, 23,6967 Appendix E Table A3. Comparison of the flow profiles obtained from experiments (with theory) actuation. The first column shows the flow profiles from the theory, assuming a sinu profile for individual pumps. The second column shows the flow profiles from the exper the actual flow response of the pumps for sinusoidal actuation. The experimental behav from the expected theory due to the threshold pressure needed for opening the valves, fl and the capacitance of the flexible tubing, which were not considered by the theory. Th umn shows both the theoretical and experimental waveforms together for a be er com the three-phase rectifier, there is half the number of ripples in the experiment (comp

Appendix C
For a one-phase rectifier pump, i = 1; ϕ 1 = 0 For a two-phase rectifier pump, i = 2; ϕ 1 = 0, ϕ 2 = π For a three-phase rectifier pump, i = 3; ϕ 1 = 0, ϕ 2 = 2π/3, ϕ 3 = 2 × 2π/3 Appendix D Figure A2. The complete experimental setup used for testing the single-and multi-phase rectifications. The microfluidic setup consists of an inlet reservoir, micropumps, a flow restrictor, a flow sensor (connected to a monitor), a microfluidic chip, and a waste reservoir. The electrical setup consists of power sources, amplifiers, RC circuits, a microcontroller, a function generator, and oscilloscopes. A function generator was used to generate a wide range of frequencies for one-phase rectification. The Arduino microcontroller was used for the low-frequency tests for single-and multiphase rectification. The computer was used to generate sine waves in Arduino software and adjust the frequency. A scissor jack and stage were used in the back-pressure experiments. A flow restrictor was used to check the performance of multi-phase rectification at very low flow rates. Figure A2. The complete experimental setup used for testing the single-and multi-phase rectifications. The microfluidic setup consists of an inlet reservoir, micropumps, a flow restrictor, a flow sensor (connected to a monitor), a microfluidic chip, and a waste reservoir. The electrical setup consists of power sources, amplifiers, RC circuits, a microcontroller, a function generator, and oscilloscopes. A function generator was used to generate a wide range of frequencies for one-phase rectification. The Arduino microcontroller was used for the low-frequency tests for single-and multi-phase rectification. The computer was used to generate sine waves in Arduino software and adjust the frequency. A scissor jack and stage were used in the back-pressure experiments. A flow restrictor was used to check the performance of multi-phase rectification at very low flow rates.

Comparison of the Flow Profiles (Theory and Experiments)
One-phase Appendix E Table A3. Comparison of the flow profiles obtained from experiments (with theory) for 0.25 Hz actuation. The first column shows the flow profiles from the theory, assuming a sinusoidal flow profile for individual pumps. The second column shows the flow profiles from the experiments with the actual flow response of the pumps for sinusoidal actuation. The experimental behavior deviates from the expected theory due to the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered by the theory. The third column shows both the theoretical and experimental waveforms together for a be er comparison. For the three-phase rectifier, there is half the number of ripples in the experiment (compared to the theory) over the same time period due to the deviation in the experimental flow profile from the theoretical values.

Comparison of the Flow Profiles (Theory and Experiments)
Two-phase Appendix E Table A3. Comparison of the flow profiles obtained from experiments (with theory) for 0.25 Hz actuation. The first column shows the flow profiles from the theory, assuming a sinusoidal flow profile for individual pumps. The second column shows the flow profiles from the experiments with the actual flow response of the pumps for sinusoidal actuation. The experimental behavior deviates from the expected theory due to the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered by the theory. The third column shows both the theoretical and experimental waveforms together for a be er comparison. For the three-phase rectifier, there is half the number of ripples in the experiment (compared to the theory) over the same time period due to the deviation in the experimental flow profile from the theoretical values.  Table A3. Comparison of the flow profiles obtained from experiments (with theory) for 0.25 Hz actuation. The first column shows the flow profiles from the theory, assuming a sinusoidal flow profile for individual pumps. The second column shows the flow profiles from the experiments with the actual flow response of the pumps for sinusoidal actuation. The experimental behavior deviates from the expected theory due to the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered by the theory. The third column shows both the theoretical and experimental waveforms together for a be er comparison. For the three-phase rectifier, there is half the number of ripples in the experiment (compared to the theory) over the same time period due to the deviation in the experimental flow profile from the theoretical values.  Table A3. Comparison of the flow profiles obtained from experiments (with theory) for 0.25 Hz actuation. The first column shows the flow profiles from the theory, assuming a sinusoidal flow profile for individual pumps. The second column shows the flow profiles from the experiments with the actual flow response of the pumps for sinusoidal actuation. The experimental behavior deviates from the expected theory due to the threshold pressure needed for opening the valves, fluid inertia, and the capacitance of the flexible tubing, which were not considered by the theory. The third column shows both the theoretical and experimental waveforms together for a be er comparison. For the three-phase rectifier, there is half the number of ripples in the experiment (compared to the theory) over the same time period due to the deviation in the experimental flow profile from the theoretical values.  Table A4. Comparison of the theoretical and experimental flow profiles and fluidic ripple factors of the rectifiers obtained for 0.5 Hz of actuation.

Waveform Graphs
Initial Waveforms in the Model Obtained Experimental (Zoomed in from Figure 6

) Waveforms
One-phase Two-phase Three-phase